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Abstract—This letter presents a new algorithm for the precise 
estimation of the frequency of a complex exponential signal in ad- 
ditive, complex, white Gaussian noise. The discrete Fourier trans- 
form (DFT)-based algorithm performs a frequency interpolation 
on the results of an N point complex fast Fourier transform. For 
large N and large signal to noise ratio, the frequency estimation 
error variance obtained is 0.063 dB above the Cramer—Rao Bound. 
The algorithm has low computational complexity and is well suited 
for real time applications. 


Index Terms—Discrete Fourier transform (DFT), fast Fourier 
transform (FFT), frequency estimation, interpolation. 


I. INTRODUCTION 


HERE has been substantial prior work in the use of the 

fast Fourier transform (FFT) to estimate the frequency of 
a time sampled complex exponential signal in additive white 
Gaussian noise [1]-[4]. This letter presents a new discrete 
Fourier transform (DFT)-based method to obtain a very accu- 
rate frequency estimate [5]. In particular, the coarse estimate 
is the frequency corresponding to the maximum amplitude 
FFT coefficient. Recursion can be done either by frequency 
translating the signal or frequency modifying DFT coefficients. 
Frequency modifying the DFT coefficients is suggested in [7]. 
A frequency error function is defined in terms of two modified 
DFT coefficients. The first interpolated frequency estimate is 
the initial frequency estimate plus the output of the frequency 
error function. Two new modified DFT coefficients are then 
obtained. The second interpolated frequency is then the first 
interpolated frequency plus the output of the frequency error 
function. This iterative computation can be continued until 
a fixed point solution is obtained in which the input and the 
output of the frequency error function are zero. Convergence 
occurs for the third frequency interpolation for the particular 
frequency error function considered in this letter. 


II. COARSE FREQUENCY ESTIMATION 
The received signal plus noise, 7[n], is given by 


r[n] = s[n] + n[n], forn=0,1,2,...,N—-1 (D 


where s[n] is a complex exponential signal with frequency f 
given by s[n] = Ae??*/"Ts and nn] is a sequence of inde- 
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pendent, identically distributed complex Gaussian random vari- 
ables with zero mean and variance o?. It is desired to process 
r[n],n = 0,1,..., N — 1 to obtain an estimate of f, where 
f is a fixed and unknown parameter, f € [0, f,]. The sampling 
frequency is f,, and T, = 1/f,. The signal-to-noise ratio (SNR) 
is defined as 


SNR = A?/o?. (2) 


Rife and Boorstyn [1] described a technique for the estimation 

frequency using the FFT. The frequency corresponding to the 

maximum amplitude FFT coefficient is chosen as a frequency 

quantized approximation to the maximum-likelihood estimate. 
Define 


r(0) y(0] 
r[1] Y([1 
t= and Y= 
r[N 1] y[N-]] 


where, Y = FFT(r) and FFT(.) is the N point complex FFT 
operator. Following Rife and Boorstyn [1], a coarse frequency 
estimate, f,, may be obtained from, 


Imax = max" [¥ [A]: 0 << N=] 


r Kimax 
fo =a Ss (3) 


Assuming that the SNR is sufficiently high, it is highly probable 
that f € [fo — (fs/2N), fo + (fs/2N)]. This is an above 
threshold condition [1]. A fine interpolation may be obtained to 
improve the frequency accuracy. 


II. DFT INTERPOLATION 


Define the modified DFT coefficients as a@ and @ as 


IN ee 
= _ 4 max 5 
an-Y (nas 5) = ; r[n] exp (- jon “meat | (4) 


n=0 
N-1 1 
kimax OY 
B-Y (Kms 5) = r[n] exp (-i2m vzect . (5) 
2 = N 


A number of possible functional forms, which map a and B 
to an estimate of the error in f,, are possible. The functional 
mappings contain a particular class [5] characterized by 


Afo(r) = ake: [67 = lal” fs 


= f 0. 6 
2Ny [a + lal tease <A 
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Fig. 1. Normalized frequency discriminant for y = 2 as a function of 


normalized frequency error. 


The parameter + defines the shape of the frequency discriminant 
function as a function of 


(Sst) 


The application of the frequency discriminant function with 
7 = 2 is investigated in this letter. This function is 


_ 1 (62 -laP 
AN |p? + lal? 7 


Fig. lshows Af, (r)/(f,/N) asafunction of (f/(fs/N)) —kmax 
for y = 2. 

From Fig. 1, itis obvious that for |(f/(fs/N)) — kmax| < 0-1, 
the discriminant yields a good approximation to the actual 
frequency error. This characteristic may be effectively utilized 
through the use of recursion. 


Afo(r) (7) 


IV. RECURSION ALGORITHM 


The algorithm is defined as 


- kimax 
fo = mers (8) 


For m = 0, 1, 2, ..., define 
N-1 a 1 
Am = e r[n] exp (~i2nn( 4 - zr) (9) 
n=0 
N-1 
cs 


m= 2 r[n] en (—i2en( jf + xx) 

D. A: 2 
Ape(e 2s leche (11) 
Fat =fm +t Afm(r)- (12) 


AN |B? + lem|? 


The interpolated frequency estimate is tsi 
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V. TRUNCATION OF THE RECURSION ALGORITHM 


foo is the implicit solution of, A fo(r) = 0. Fae is the 
m + 1th interpolation. In practice, te is sufficiently close to 
foo to be used as the interpolated frequency estimate. When 
\(f/(fs/N)) — kmax| > 0.25, f, will have reasonably large 
error. However, |(f — f1)/(fs/N)| < 0.25 and therefore fy 
will have small error and A f will be close to zero. Therefore, 
fs will be essentially the implicit solution to Af, (r) = 0. The 
recursion may be truncated to the evaluation of fas 


VI. REDUCTION OF THE RMS ERROR USING RECURSION 


The rms frequency estimation error due to additive noise, 
using the proposed frequency discriminant, decreases sharply 
for (f — - )/(fs/N) & 0. Since the iterative solution adaptively 
produces a frequency estimate, which is close to this condition, 
the recursive use of the discriminant produces a frequency esti- 
mate with small rms error. By means of the recursion, the dis- 
criminant converges to a minimum rms error condition. 


VII. PERFORMANCE OF THE ALGORITHM 
Using the Taylor series [6], at high signal to noise ratios, the 
performance of the algorithm has been shown to be [5], 
N sin? (4) tan? (34) 
4(SNR) 72 


oe (13) 


where, oF is the mean square value of ( f= f)T; and f isa 
uniform random variable in the interval [0, fs]. 

The Cramer—Rao lower bound (CRLB) on the unbiased fre- 
quency estimator mean square error [1] is 


6 


2 = 
CCRLB ~ (97)2N(N2 — 1)(SNR)’ me 


Therefore, the performance of the frequency estimation algo- 
rithm compared to the CRLB is, 
2 2/72 +59 2 
of _ N?(N? —1)sin (sy) tan? (47) (15) 
2 _ 6 ‘ 
OCRLB 


For large SNR and large N 


oF a4 
li = 101 — |= 0.0633 dB 
> g ORLB oe ( 96 ) 
N=o (16) 


This limit was also obtained by a different method in [7]. 


VIII. SIMULATION OF THE PERFORMANCE 


The performance of the frequency estimator was obtained by 
simulation. 

Figs. 2 and 3 show the rms estimation error, of, as a function 
of the SNR in decibels for the cases of a 16-point and a 1024- 
point FFT, respectively. The results agree with the analysis. 
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Fig. 2. The rms frequency estimation error as a function of the SNR in 
decibels. Results are shown of one, two, and three iterations. The results are 
compared to the CRLB. The FFT length is V = 16. 


IX. CONCLUSION 


An algorithm for frequency estimation has been introduced. 
The high SNR performance of the algorithm is 0.063 dB above 
the CRLB. The algorithm has low computational complexity 
and is suitable for real time digital signal processing applica- 
tions including communications, radar and sonar. 
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